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ABSTRACT 


Two luminaires were evaluated to determine the light flux density pattern 
on a horizontal plane surface. NASA supplied both luminaires; one was made by 
NASA and the other is commercially available. Tests were made for three 
combinations of luminaire height and luminaire lens material using the NASA 
luminaire; only one configuration of the commercial luminaire was tested. 
Measurements were made using four sensors with different wavelength range 

capabilities. The data are presented in graphical, tabular, and magnetic disk 
formats . 


RESEARCH FACILITIES 

The Luminaire Evaluation Chamber (LECH) in the Agricultural Engineering 
Building at Penn State University was used to collect data. The chamber 
contains a stand to mount the luminaire at the desired height and to rotate 
the luminaire about a vertical axis. A linear track on the floor of the 
chamber is used to position the light sensor at desired radial distances from 
the axis of rotation of the luminaire. The rotation of the luminaire and the 
movement of the sensor are effected by stepping motors controlled by a 
microcomputer, fhe microcomputer simultaneously collects data from the light 
sensor and records the data on a magnetic disk. The luminaire stand and 
sensor track are enclosed in a structure 7.7 x 2.5 x 5.0 m high, constructed 
with waferboard interior surfaces, including the floor, which have been coated 
with flat black latex paint. 

The sensor track runs the long direction of the structure and the 
luminaire stand is located at one end. The chamber also contains an 
analog-to-digital device for conditioning signals to and from the computer, 
two stepping motor controllers, a sensor signal amplifier, a voltage regulator 
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to control voltage to the luminaire, two box fang to circulate the air, and a 
thermocouple air temperature recorder. A black cotton curtain wag used to 
shield the sensor from any light produced by the computer or instrumentation. 
The layout of the structure is shown in Figure 1. The chamber is illuminated 
with four 40-W fluorescent tubes mounted to the ceiling. These lamps were off 
during data collection. 


SCOPE OF TESTS 

A NASA designed luminaire and a commercially available lumimaire, both 
provided by NASA, were evaluated in the Luminaire Evaluation Chamber. The 
NASA lumimaire was evaluated with a glass lens and also with a plastic 
(polymethylmethacrylate; Plexiglass) lens. The commercial luminaire had no 
lens* (The te ™ lens is used here to mean a transparent bottom cover.) All 
tests were done at a luminaire mounting height of 0.5 m above the sensors, 
except for the NASA luminaire with the glass lens which was also tested at a 
height of 1.0 m. Considering each combination of luminaire, lens material, 
and height as a "luminaire configuration", a total of four luminaire 
configurations were tested. 

For each luminaire configuration, light measurements were made with four 
sensors. 1) a Li-Cor LI-190SB (Model: QUANTUM, S/N: Q6724, Calibration date: 
89 Aug 11} which measures photosynthetic photon flux density Integrated over 
the 400-700 nm wavelength range. 2) an Eppley Precision Spectral Pyranometer 
(Model: PSP, S/N: 26456F3, Calibration date: 89 Sept 18) with a Schott RG805 
glass hemisphere which measures energy flux density integrated over the 
800-3000 nm range. A 3 mm long scratch was noted near the top of the glass 
hemisphere of this sensor. 3) an Eppley Precision Infrared Radiometer (Model: 
PIR, S/N: 26407F3, Calibration date: 89 Sept 12) which 


measures energy flux 


density integrated over the 3,000-50,000 nm range. And 4) a Li-Cor Portable 
Spectroradiometer (Model: LI- 1 800, S/N: PRS-159, Calibration date: 89 Jan) 
which measures energy flux density over 2 nm wavelength increments in the 
300-850 nm range. These four sensors will be referred to as the Quantum, PSP, 
PIR, and PRS sensors, r, 3 pectively. 

SENSOR MEASUREMENT GRID 

Measurements for all sensors were made within what represents a 2.5 m 
radius circle on a horizontal plane beneath the luminaire. Measurement 
locations were reached by moving the sensor incrementally along the sensor 
track and sequentially rotating the luminaire. The Quantum, PSP, and PIR 
sensors were mounted to a computer-controlled cart and measurement locations 
made a polar grid with radial increments of 0.05 m and angular increments of 
3.6 degrees. The PRS sensor was too heavy to be moved by the cart and was, 
therefore, moved by hand. PRS measurements were made at radial increments of 
0.5 m and angular increments of 45 degrees. 

The center of the polar grid was directly beneath the center of rotation 
of the luminaire. The center of rotation for the NASA luminaire was along the 
axis of symmetry, 597 mm from the surface behind the air filter. This 
location was experimentally chosen so the axis of rotation would pass through 
the approximate center of gravity of the luminaire. Thus, the forces on the 
luminaire stand, and therefore stand deflection, remained constant at all 
angles of rotation. 

The NASA luminaire was mounted onto the stand which was then adjusted so 
its axis of rotation was vertical; next, the luminaire mounting was adjusted 
so the luminaire opening was horizontal. A 0.5 m long level laid on top of 
the glass lens was used to observe that the opening remained horizontal at all 


angles of rotation. The lamp ballasts and related equipment were mounted 


remote from the luminaire. 

The same mounting procedure was used for the commercial luminaire except 
the axis of rotation was chosen to pass through the center of the luminaire 
opening and the luminaire was adjusted to horizontal by holding the level 
across the bottom of the fixture. Also, the ballast was left attached to the 
luminaire. 

The zero angle for the polar grid is along a line parallel to the nominal 
lamp axis, with a ray extending from the base to the tip of the lamp being in 
a positive direction. The positive angular direction is counterclockwise when 
viewed from the top. This means, for example, that the luminaire was rotated 
clockwise 3.6 degrees in order to take the readings along the +3.6 degree r ay. 
The order in which measurements were taken using the computer controlled cart 
was determined by the program COLLECT. BAS. (Appendix B contains all 
user-generated programs.) The order is Illustrated in Figure 2. This 
measuring scheme gave redundant measurements at the origin and points along 
the first ray which was both the zero and 360 degree ray. Multiple readings 
at a single point may be compared to determine the stability of the 
measurements . 

The order of measurements with the PRS sensor did not follow a rigid 
pattern except one radius was chosen, then measurements were made at all 
angles for that radius before another radius was selected. At the first 
selected radius the order of data collection was from zero to 360 degrees; at 
the second selected radius the order was 360 to zero; etc. The luminaire was 
rotated using the computer programs CCW^5.BAS and CW^S.BAS, respectively, for 
each radius. The data files for the PRS sensor contain the time of day so the 
exact order of data collection can be reconstructed. 
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DATA COLLECTION PROCEDURE 

For ail tests the luminaire was allowed to staDilize for a minimum of two 
hours after turning it on. Likewise, the amplifier was allowed to warm-up for 
a minimum of two hours before using it. The computer-controlled cart was 
moved to a radial position of about 1.25 meters. Before each test the 
appropriate Quantum, PSP, or PIR sensor was attached to the cart and leveled. 
(The Quantum sensor had no leveling bubble and, therefore, was not precisely 
leveled . ) 

The analog to digital conversion software was initiated and then the 
program COLLECT. BAS was run with the stepping motors switched off. The sensor 
signal cable was disconnected at the sensor, a 100 ohm resistor was installed 
across the two signal conductors, and the amplifier was adjusted so the 
readings taken by the COLLECT. BAS program (and displayed on the computer 
screen) were approximately zero. The computer program was then stopped. 

The cable was reconnected to the sensor, the cart was manually moved to 
its zero position, the luminaire was manually moved to its zero angle 
position, the stepping motors were switched on, the temperature recorder was 
started to collect air temperature every 15 minutes. COLLECT. BAS was 
restarted to begin data collection; the program allowed a 20 second delay 
before the first reading and in that time the operator left the chamber. 

DATA ON MAGNETIC DISKS 

The data is contained on 1H floppy disks (5-1A, ASC T I, MS-DOS). A 
fifteenth disk contains computer programs used in this study. Two of the 
fourteen disks contain data from the PRS sensor and are labeled PRS1 and PRS2. 
The files on these two disks are described in the DATA IN TABULAR FORMAT 


section. The other twelve disks contain data from the Quantum, PSP, and PIR 
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sensors . They are labeled with an seven character string beginning with the 
letters Q, S, or I to indicate which sensor. The second and third characters 
in the string indicate the height of the luminaire above the sensor: 05 is 

0.5 m and 10 is 1.0 m. The fourth character represents the luminaire 
configuration: NASA luminaire with Glass lens, NASA luminaire with Plastic 

lens, or Commercial luminaire. The fifth, sixth, and seventh characters 
represent the month and day the data were collected: N28 is November 28. 

Each of the twelve disks contain three data files labeled with an eight 
character file name with a three character extension. The first seven 
characters in the file name are as described immediately above. The eighth 
character is 0 if the file contains sensor readings from the COLLECT. BAS 
program. The eighth character is C if the sensor reading has been Changed or 
modified using the program LUMCHu.FOR as described below. The extension .RAW 
indicates the file was produced by the program COLLECT. BAS while the extension 
.DAT indicates it was produced by the program LUMDAT.FOR. The extension .GRD 
indicates the file was produced by the Golden Software Inc. Surfer software 
program GRID.EXE. 

The file on each disk with ;he extension .RAW contains the data as taken 
from the Luminaire Evaluation Chamber. Each line of data contains two numbers. 
The first is the order number of the sensor reading (see Figure 2) and the 
second is the sensor reading, converted to the appropriate units. The program 
COLLECT. BAS converts the millivolt output of the sensor to umol/m 2 s for the 
Quantum sensor and W/m 2 for the PSP and PIR sensors. 

The file on each disk with the extension .DAT contains data from the .RAW 
file on the same disk, but which has been modified using programs LUMCHG.FOR 
and LUMDAT.FOR. LUMCHG.FOR looks at the data from each ray in the polar grid 
and subtracts the sensor reading at the 2.5 m radius from all the other 
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readings along that ray. The purpose was primarily to differentiate between 
the effects of the luminaire and the chamber for the PIR data. The secondary 
reason was to minimize the effect of signal drift in the instrumentation, and 
for that reason the program was used on the Quantum and PSP data as well. 
LUMCHG.FOR outputs to an intermediate file which LUMDAT. FOR reads. 

LUMDAT.FOR assigns cartesian coordinates to the data based on their order in 
the file and outputs to a .DAT file. Each line of a .DAT file contains an 
X-coordinate , Y-coordlnate, and the 3 ensor reading. The coordinates ha ” 3 
units of meters and the sensor reading has the same units as before. 

Golden Software Inc. SURFER software program GRID.EXE takes the data from 
the polar grid (albeit with cartesian coordinates) and interplates and 
extrapolates the data to produce data in a uniform cartesian coordinates grid. 
The new grid extends from -2.5 to 2.5 meters in both directions with nodes 
every 0.05 m for a total grid size of 101 x 101 nodes. GRID.EXE contains many 
options for data conversion and the options used in this study are specified 
by the file GRIDNASA.CMD . GRID.EXE outputs the new grid to a .GRD file. The 
ten lines of information at the top of the .GRD file are plotting Information 
for the other SURFER programs SURF.EXE and T0P0.EXE. The remaining lines are 
the Interpolated data arranged serially for node locations In the grid read 
left to right, line by line, starting at the node (-50,-50) which Is 
equivalent to location (- 2.5 m - 2.5 m). 

DATA IN TABULAR FORMAT 

The data for each of the four luminaire conf ' gurations are together in a 
separate volume. Each volume contains the data for all four sensors. The 
first three sections of each volume contain the data from the Quantum, PSP, 
and PIR sensors and the data In these three sections are arranged In the same 
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format. On the title page of each section is the chamber air temperature 
range during data collection. The first two columns of data in each section 
give the position of the sensor reading in polar coordinates, and the second 
two columns are the cartesian equivalents of those coordinates. The last 
column is the sensor reading. These tables were printed using the program 
LUMPRI.FOR operating on files produced by the program COLLECT. BAS (l.e., they 
were not modified using the program LUMCHG.FOR ) . 

The fourth section contains data from the PRS1 and/or PRS2 floppy disks, 
collected with the PRS sensor. This section contains 54 two-page files, one 
for each sensor position. The sensor position is coded into the three 
character file name. The first character indicates the luminaire 
configuration; NASA luminaire with Glass lens at 0.5 m height, NASA luminaire 
with Glass lens at 1.0 m Height, NASA luminaire with plastic 
(Polymethylmethacrylate ) lens, and Conventional luminaire. The second 
character in the file name is a number which represents the radius in the 
polar grid. This number multiplied by 0.5 gives the radius in meters. The 
third character is a number which represents the angle in the polar grid. 

This number multiplied by 45 gives the angle in degrees. 

The first seven lines of each file gives summary Information which is 
self explanatory. The rema'-.lng lines are a serial listing of the data with 
the wavelength followed by the sensor reading at that wavelength. Wavelength 
units are nm and sensor reading units are micromoles per square meter second. 
These files began as files created by the PRS instrument. The PRS files were 
transformed to quantum units using the Li-Cor software program LI1300 
following the FILE OPERATIONS, VIEW, and QUANTUM TRANSFORM options and 
suboptions. The file thus created was printed using the LIPRI.FOR program. 
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DATA IN GRAPHICAL FORMAT 

Data from the Quantum, PSP, and PIR sensors for each of the four 
luminaire configurations are presented graphically in Appendix A. Two types 
of graphs were drawn for each data set: orthographic and contour. Both type* 
of graphs were created using the Golden Software Inc. SURFER programs 
operating on files with .GRD extensions. The origin of the .GRD files is 
described in the section DATA ON MAGNETIC DISKS. The orthographic graphs were 
produced using the program SURF.EXE in conjunction with the command file 
SURFNASA.CMD , while the contour graphs were produced using the program 
T0P0.EXE in conjunction with the command file T0P0NASA.CMD. and the file 
NASA3.LVL. 
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COMPUTER PROGRAMS 

Numerous "user-generated" computer programs were written or modified 
specifically for this study and copies are provided herein. Several 
commercial programs were also used as noted throughout this report, but copies 
of these programs are not listed here. The user-generated programs are 
provided both on a magnetic disk labeled "Computer Programs" and as hard copy 
listings in Appendix B. These programs can be categorized as one of three 
types: GWBASIC, WATF0R77 FORTRAN, and Golden Software. 

The GWBASIC programs are those used to operate the Luminaire Evaluation 
Chamber. These programs are COLLECT. BAS, CCW45.8AS, and CW45 . BAS . 

COLLECT. BAS was modified during the study. When used for the NASA luminaire, 
all sensor readings less than zero were set to zero. Later, when used ro r the 
commercial luminaire, negative values were allowed to remain. The effects of 
this change are believed to be Insignificant. 
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GRIDNASA.CMD , SURFNASA.CMD , AND TOPONASA.CMD are not true programs; they 
are command fiiC3 specifying some of the options selected within the Golden 
Software Inc. SURFER programs. These files were generated by and may be used 
by the appropriate SURFER programs. NASA3.LVL specifies which contour lines 
will be graphed when using TOPONASA.CMD and SET21 .STM is the symbol set used 
to label the graphs in Appendix A. It contains the micron and other special 
symbols. Unlike the other files, SET21 .SYM is provided only on magnetic disk. 

Programs written in WATFOR77 FORTRAN were used to manipulate and print 
data flies. These programs Include LUMCHG. FOR , LUMDAT . FOR , LUMP RI. FOR, and 
LIPRI.FOR. The OPEN commands in these programs were edited to read and write 
to specific data flies. The program PLANCK. FOR is also included; it is used 
calculate the flux density from a blackbody over a specified wavelength range 
based on Planck's law. This program was not used directly in this study, but 
ratner , was used to check the output of the PIR sensor. An executable copy of 
this program, PLANCK.EXE, Is provided on magnetic disk. 
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Figure A13. NASA luminaire, 1.0 m height, glass lens. Quantum sensor , contour 
line units are mol/m 2 s. 









Figure A16, NASA luminaire, 0.5 m height, glass lens, Quantum sensor, contour 
contour line units are pmol/m 2 s. 
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Figure A2^. Commercial luminaire, 0.5 m hetght, PIR sensor, contour 
lines units are W/m^. 





Appendix B. COMPUTER PROGRAMS 
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1 ’ ********************************************** ********** **** 
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3 ’ CCW45.BAS 

4 ’ ************************************************************ 

5 ’ ************************************************************ 
10 DATA &hb8 , &h59 ,&h47 ,&hcd,&h60 ,&h90 ,&h90 , &hca,&h06 , &h0 , 0 

20 DEF SEG 

30 APROG 1$= SPACES ( 15) 

40 APROG 2 S= SPACES (15) 

50 A%=VARPTR(APR0G1$) : AM1=PEEK( A%+1 )+PEEK( A%+2 ) *256 
60 RESTORE 10 
70 FOR A=AM1 TO AM1+10 
80 READ A% : POKE A,A% 

90 NEXT A 

100 A%=VARPTR( APROG2S) : AM2=PEEK(A%+1 )+PEEK( A%+2 )*256 

110 RESTORE 10 

120 FOR A=AM2 TO AM2+10 

130 READ A%: POKE A,A% 

140 NEXT A 

150 POKE AM2+8 , 16 

1000 C$ = "S" + CHR$(0) :A%(2)=l:A%(3)=i : CALL AMI ( AX(0 ) , B(0) , C$) 
1500 FOR I = 1 TO 1250 

2000 C$ = "0" + CHR$ ( 0 ) : A% ( 2 ) = 1 : CALL AMI ( A%(0 ) , B( 0 ) , C$) 

2010 C$ = "0" + CHRS ( 0 ) : A% ( 2 ) =0 : CALL AMI ( A%( 0 ) , B(0 ) , C$ ) 

2050 NEXT I 


11 > *********************************************************** 

2 * x*********************************************************** 

3 ’ CW45.BAS 

4 9 ************************************************************ 

5 » ************************************************************ 

10 DATA &hb8 , &h5 9 , &h4 7 , &hcd , &h60 , &H90 , &h90 , &hca , &h0 6 , &h0 , 0 
20 DEF SEG 

30 APR0G1$=SPACE$( 15) 

40 APROG2$=SPACE$ ( 15 ) 

50 A%= VARPTR ( APROG 1 $ ) : AM1=PEEK( A%+1 )+PEEK( A%+2)*256 
60 RESTORE 10 
70 FOR A=AM1 TO AM1+10 
80 READ A% : POKE A, A% 

90 NEXT A 

100 A%=VARPTR(APROG2$) : AM2=PEEK(A%+l)+PEEK(A%+2)*256 

110 RESTORE 10 

120 FOR A=AM2 TO AM2+10 

130 READ A%: POKE A,A% 

140 NEXT A 

150 POKE AM2+8.16 

1000 C$ = "S“ + CHR$(0) : A%(2)=1: A%(3)=1=CALL AM1(A%(0) ,B(0) ,C$) 
1500 FOR I = 1 TO 1250 

2000 C$ = "O'* + CHR$ ( 0 ) : A% ( 3 ) = 1 : CALL AM1(A%(0) ,B(0) ,C$) 

2010 C$ = ”0" + CHR$ ( 0 ) : A% ( 3 ) =0 : CALL AMI ( A% ( 0 ) , B ( 0 ) , C$ ) 

2050 NEXT I 




col =10 

gmampow =2 

gmamcen =2 

gmamnx =2 

co 2=112 

gmsmny =2 

co3=15 

gmbinps" “ 

snow =7 

gmbrl=l 

bios =n 

gmbr2=32767 

path=" " 

gmbcl=l 

gm=i 

gmbc2=32767 

grmpower=n 

gmbout = ” OOT . GRD 

gnunaxerrora=n 

gmboutfmf =b 

gtnmaxi tera =a 

gmbblks " " 

am=n 

gmminpls” ” 

maxrot =10 

gmminp 2 =” ” 

maxrad=l. 5 

gnunout = " OUT . GRD ' 

gronamo= “ M 

gmmequs " " 

groutfmt=a 

gmmlrl=l 

gridupflag=i 

gmmlr2=32767 

grddistx=0.05 

gmmlclsl 

grddiaty=0.05 

gmmlc2=32767 

gr<. 7 oqx=l 

gmm 2 rl=l 

gxmin=- 2 . 5 

gnm2r2=32767 

gxmaxs 2 . 5 

gmm 2 cl=l 

gymin=- 2 . 5 

gmm2c2=32767 

gymax= 2 . 5 

gmmoba-b 

gldiac=n 

xformequ="D = A 

grxcol=l 

gwkasf=" ” 

grycol =2 

grinames ” " 

gr 2 col =3 

grinfmt=d 

gfunequ ='*2 = x*x + y*y" 

maxcol=3 


gf unxmax=25 

gf unxinc=l 

gf unyrain=0 

gf unymax=25 

gfunyinc=l 

gf unog="OUT . GRD" 

gfunab=b 

gmainps '* " 

graaouts "OUT . GRD " 

gmarl=l 

gmar2=32767 

gmscl=l 

gmac2=32767 

gmabaout=b 

gmanx=2 

gmsny=2 

gmamtype=a 




1*7 

♦♦at*************:*#****:***************#****#**********:*******:***** 

HX***********************************-*#************************** 

SURFNASA.CMD 

a***************************************************************: * 
****************#**********#*****:M»****************************** 


col =10 

co 2=112 

co3=15 

acrtype=a 

acrwid=14 

acrpoas-2, -0 . 75 

snow =7 

bioa =n 

egapals" " 
paths" " 
plotunit=i 
innames" " 

inrl=l 

incl=l 

inr2=32767 

inc2=32767 

pro j=o 

itb=u 

ihesy 

angh=225 

angv=30 

oeyediatsa 

diratr="XY" 

plotpen=l 

cmin=-l 

cmax=ll 

cint=l 

levnames " " 

drawborder=y 

czafilos 

ipb=y 

ipbl=n 

basepen=l 

obaaeht =0 

titatr=" " 

titaymnumsO 

tpoa=a 

titangsO 

tithtsO . 175 

titpensl 

legatatsa 

lpoa=a 

legpen=l 

textf ilemode=a 

textf ile=" 

f ontO = " SET2 1 . SYM " 


f ontl="DEFAULT . SYM" 
font2s”DEF AULT. SYM" 
font 3= "DEFAULT . SYM" 
font4= "DEFAULT. SYM" 
axiatitleaym =0 
axiatitlepen = 1 
axistitlesiz=0 . 15 
axistitleap= 0 . 2 
xplotax-y 
xaxispen=l 
xaxsaym =0 
xoplane=a 
xoticdiat=0 . 25 
xotlfreq=4 
xtlfmt=f 
xt lndig = 1 
xtlaize= 0 . 1 
xtlangle =0 
xtlspsO . 1 

xtitles"LAMP AXIS (m)" 

xtitaiz=0 . 15 

xtitapsO . 2 

yplotaxsy 

yaxispen=l 

yaxssymsO 

yoplane=a 

yoticdiat=0 . 25 

yotlfreq=4 

ytlfmt=f 

ytlndigsl 

ytlaize= 0 . 1 

ytlangle =0 

ytlapsO . 1 

ytitlas" TRANSVERSE AXIS (m) 
ytitsiz=0 . 15 
ytitapsO . 2 
zplotax=y 

zaxispensl 

zaxaaymsO 

zoplane=a 

zoticdiatsl 

zotlf req=5 

ztlfmtsf 

ztlndig =0 

ztlaizesO . 1 

ztlanglesO 
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z"tlap=0 . 1 
2w) " 

ztitaiz=0 . 15 
ztitap=0 . 2 
zatart=0 
zend= 10 
aidelen=5 
map7eq>.= l 
ozacale=a 
x 7 namo=" " 

X7pen=l 

poatname="A: I10GN070.DAT" 

postcola='' 1, 2, 3" 

poat£mt=c 

poatndig=2 

poat33nnnum=0 

plht=0 . 1 

plang=0 

labline=l 

plpen=l 

plpoa=a 

poatf ilemode=a 
poatf ile=" ” 

outname="I10GN070 .PLT” 

outacale=l 

pagepoa=l, 1 

pftype=b 

outmode=o 

outndig=3 

callplot =7 

plotapec= "PLOT . EXE" 


X**************************************************************** 

***************************************************************** 

TOPONASA.CMD 

***************************************************************** 

x**************************************************************** 


col =10 

co 2=112 

co3=15 

scrtype=a 

scrwid=14 

scrpos=-l ,-0.75 

snow =7 

bios =n 

egapal=" " 

path= ” ” 

plotunit=i 

inname =" " 

inrl=l 

incl=l 

inr2=32767 

inc2=32767 

cmin =10 

cmax=500 

cint= 

1 evname = " NAS A3 . LVL" 

maplen=5 . 5 

mapy eqx= 1 

titstr=" 

titaymnum =0 

titpos=a 

titang =0 

titht=0 . 175 

titpen=l 

legstat=n 

legsymnum =0 

legunits="data units" 

legpos=a 

leght=0 . 15 

legpen=l 

botstr="LAMP AXIS (m)“ 

If tstr= " TRANSVERSE AXIS (m)" 

rhtstrs" H 

topatrs" " 

btitsym =0 

btitht=0 . 175 

btitpen=l 

textf ilemode=a 

textf ile=" " 

f ontO="SET21 . SYM" 

fontls" DEFAULT. SYM" 

font 2s" DEFAULT. SYM" 


f ont 3 = " DEFAULT . SYM" 

font 4= "DEFAULT . SYM" 

borstatsy 

borticpen=l 

ticside="LRTB" 

xxtic=0 . 25 

xltf req=4 

xborangsO 

yytic=0 . 25 

yltf req=4 

yborangsO 

bcrlabfmt=f 

borndigsl 

borlabpen=l 

borht= 0 . 1 

borsymnum =0 

xyname= " " 

xypen = 1 

grdpen- 1 

grddash =0 

grdxx=n 

grdyy=n 

postnames "\LECH\I 10N0 100 . DAT" 

postcols=" 1 , 2 , 3 , 0 , 0 " 

postincl=n 

cencode=41 

pcensymnum =0 

cenhts= 0 . 1 

oenhtmin= 0 . 1 

cenhtmax= 0 . 3 

czmins=g 

czmaxssg 

cenhtzcol=3 

cenrotmode=d 

cenangsO 

cenpensl 

postfmt=c 

postndigs 2 

postsymnum =0 

plht= 0 . 1 

plangsO 

plpen=l 

plpos=a 

postf ileoode=a 
postfile=" " 
lconf req=l 
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lconpen=l 

lcondash=0 

lconhatch=0 

lconbold=0 

ulconpen=l 

ulcondaah=0 

ulconhatch=0 

ulconboid-0 

clab\mt=f 

clabndig=0 

clabpen=l 

clabht=0 . 1 

crvtol=4 

llmin=4 

lemin=0 . 2 

clabsjnnnum=0 

am3tat=n 

tenaion=2 

czaf 5 le=" “ 

outnaines " TEMP . PLT " 

out seal e=l 

pagepoa=2, 1 . 35 

pf type=b 

outmodo=o 

outndig=3 

callplot=y 

plotapec= " PLOT . EXE " 


**********************x*****x****xx****x*xx***xxx** M *xx M xx*xxxxxx 

x**************************************************************** 

NASA3 . L7L 

******x***x********x*****x***x*xxxx**x**xx*****xxxxxxxxxxxxxm M „x 

x**************************************************************^ 
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C****************X****X*******X**#**X*******X***X**xx*:fc***** 


C********x****x*****x*******x******^x*******x* 

DIMENSION N(101), R(101), C(101) 

OPEN ( UNIT=1 , r 'TLE= • a: I05CN290 . RAW' ) 
OPEN (UNIT=2, FILE= ' \LECH\I05CN29C RAW’) 
DO 400 1=1,50 
DO 100 J=1 , 51 
READ ( 1 , * ) N ( J ) ,R( J) 

DO 200 J=1 , 51 
C( J) =R( J ) -R( 51 ) 

WRITE ( 2, 500)N( J) , C( J) 

DO 300 J=1 , 51 
READ ( 1 , * ) N ( J ) ,R( J ) 

DO 400 Jsi, 51 
C( J ) =R( J ) -R( 2 ) 

WRITE ( 2 , 500 )N( J ) , C( J ) 

FORMAT ( 14 , IX, F7 . 3 ) 

END 


100 


200 


300 


400 

500 


*X****X******* 

.»c*x***x*****x* 
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£************* + ************************+********#****» *:*##** 
( 3 ******)*^:^** *************■***********+************** ********* 
C LUMDAT . FOR 

C.*************************u ******************************** 

C*********************************************************** 

open (unit=l, £ile= * \LECH\I05CN29C . RAW ) 
open (unit=2, files ’ A: I05CN29C . DAT ’ ) 

PI = 3.14159 
ANGLE =0 . 0 
H=0 . 0 

DO 400 1=1,50 
DO 200 J=1 , 51 
EXECUTE CONVERT 
200 3=H+0.05 

ANGLE = ANGLE + 2.*PI/100. 

DO 300 K=1 , 51 
H=H-0 . 05 
EXECUTE CONVERT 
300 CONTINUE 

400 ANGLE = ANGLE + 2.*PI/100. 

STOP 

REMOTE BLOCK CONVERT 
X=H*COS ( ANGLE ) 

Y=H*SIN( ANGLE) 

READ( 1 , * )SN, R 
WRITE( 2 , 900 ) X, Y.R 
900 FORMAT (F6. 3, 1X.F6.3, 1X.F7.3) 

ENDBLOCK 

END 


C LUMPRI . FOR 

C*********************************************************** 

C****************************************************** 

OPEN (UNITsl, FILE=’ A-*P. RAW* ) 

OPEN (UNIT=7 , FILE=’ (C)COMl’ ) 

PI = 3.14159 
ANGLE sO.O 
DANGLEsO . 0 
HsO.O 


WRITE( 7 , 906 ) 

WRITE( 7 , 910 ) 

DO 400 Isl,51 
DO 200 J=1,5X 
EXECUTE CONVERT 
200 HsH+0.05 

ANGLE s ANGLE + 2. *PI/100 
DANGLEs DANGLE +3.8 
WRITE( 7,905) 

WRITE(7, 906) 

WRITE( 7 , 910 ) 

DO 300 Ksl,5l 
HsH-0.05 
EXECUTE CONVERT 
300 CONTINUE 


ANGLE s ANGLE + 2. *PI/100 
DANGLE = DANGLE +3.6 
WRITE( 7 , 905 ) 

WRITE(7, 906) 

400 WRITE( 7,910) 

STOP 


REMOTE BLOCK CONVERT 
X=H*COS (ANGLE) 

YsH*SIN( ANGLE) 

READ ( 1 , * ) SN . R 

WRITE( 7 , 900 ) DANGLE, H.X.Y.R 

VM'HPT nrir 

900 FORMAT ( ’ M7X.F12. 1.F12.2.3F12.3) 

905 FORMAT ( ’ 1’ ) 

906 FORMAT ( ’ - ’ ) 

910 FORMAT ( ' ’,217, 'ANGLE (deg) 2X, ’RADIO: 
1,3X, ’FLUX DENSITY’ ) 

END 


(m) ’ ,4X, ’X (a) ’ ,7X, ’Y 


(m 


C*******^**j(.******»****#*** # * ¥ ,****** # ** ####J)tJ({3(tJ(t<tj({3(tj||#j(tj)c 
C LIPRI . FOR 


C*******************************************#*****^*^***** 
G************™*****************#*********####*#********** 
CHARACTER LINE*24, FN*20 
DIMENSION N(5) , V(5) 

PRINT *, 'FILE NAME (EG., A: MOO . PRN) : ’ 

READ (6,800) FN 

C OPEN (UNIT=1, FILE=’\LI1800\MOO.PRN' ) 

OPEN (UNIT=1 , FILE=FN) 

OPEN (UNIT=7, FILE=’ (C)COMl’ ) 

WRITE (7, 905) 

DO 100 1=1,7 


READ (1,800) LINE 
100 WRITE(7,850)LINE 
DO 300 J= 1 , 40 
DO 200 K=1 , 5 
200 READ(1,*) N(K) , V(K) 

300 WRITE( 7 , 900 ) (N(K) , V(K) ,K=1, 5) 
WRITE(7, 910) 

WRITE(7, 905) 

DO 500 J=l, 15 


DO 400 K=1 , 5 
400 READ(1,*) N(K) , V(K) 

500 WRITE( 7 , 900 ) (N(K) , V(K) ,K=1, 5) 
READU,*) N(1),V(1) 

WRITE( 7 , 900 ) N( 1 ) , V( 1 ) 

WRITE(7, 910) 

STOP 

800 FORMAT (A) 

850 FORMAT ( ’ ' ,11X,A24) 

900 FORMAT ( ’ ' , 10X, 5(14, E10 . 3) ) 

905 FORMAT ( ) 

910 FORMAT ( ’ 1' ) 

END 
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C*******************************************#********^**^* 

C*****jjc*****«j******* ## * #<# * ## ^ ##<# ^ #<w####jfcj(cj)cj(tj){###jfc 

U PLANCK. FOR 

%*********************************************************** 

c***************************************************^**^*** 

REAL K 
C=3 . 0E08 

C SPEED OF LIGHT IS 3.0E08 M/S 

PI=3 . 1416 
H=6 . 62E-34 
K=1 . 30E-23 

n^rwm * ’ ’THIS PROGRAM FINDS RADIANT ENERGY FLUX DENSITY 
PRINT *' BLACKBODY INTEGRATED OVER A WAVELENGTH RANGE’ 
FRINT * , * 9 

PRINT *, 'INPUT BLACKBODY TEMPERATURE (C)’ 

READ (5,*) T 
T=T+273 . 


PRINT *, ’ ’ 

BBI NT INPUT LOWER RANGE OF WAVELENGTH (NM)’ 

W1=W1 *' 1 . 0E-09 

F1=C/W1 

F=F1 


WF = 2 . 0 *P I *H/C** 2 *F** 3 / ( EXP ( H*F/K/T ) - 1 0 ) 

PRINT * t ’ EQUIVALENT FREQUENCY (CYCLES/SEC) = ’, FI. WF 
ruINT ^ ) 9 

PRINT *, ’INPUT UPPER RANGE WAVELENGTHS (NM)’ 

READ (5,*) W2 
W2=W2 * 1.0E-09 
F2=C/W2 
F=F2 


FROM A ’ 


WF=2 . 0*PI *H/C**2*F**3/ ( EXP ( H*F/K/T ) - 1 0 ) 

PRINT *, ’EQUIVALENT FREQUENCY (CYCLES/SEC) =’, F2, WF 
F HINT ^ ) 9 

READ T (5 , *) I N P0T NDMBER ° F INTERVALS FOR NUMERICAL INTEGRATION’ 
W=0 . 0 


WF1=2 . 0*PI*H/C**2*F1**3/(EXP(H*F1/K/T) -1 0 
KwF1/2*!5f ^ 2, ' tF2 * 3k3/(EX P (H * F2 / K /T)-1.0 
do iooo’ns=i,n-i 
F=F2+(F1-F2)*NS/N 


\ 

) 

) 


1000 w- W +WF* PI * H/C **2*F**3 / ( EXP ( H*F/K/T ) - 1 . 0 ) 
W=W*(F1-F2)/N 

PRINT *, ’TOTAL FLUX DENSITY = \W, ’ WATTS/M2’ 


